The chiral and angular momentum content 

of the p-meson* 

L. Ya. Glozman, C. B. Lang, and M. Limmer 

Institut fiir Physik, FB Theoretische Physik, 
Universitat Graz, A-8010 Graz, Austria 

September 16, 2009 



Abstract 

It is possible to define and calculate in a gauge-invariant manner 
the chiral as well as the partial wave content of the quark-antiquark 
Fock component of a meson in the infrared, where mass is generated. 
Using the variational method and a set of interpolators that span a 
complete chiral basis we extract in a lattice QCD Monte Carlo sim- 
ulation with rij = 2 dynamical light quarks the orbital angular mo- 
mentum and spin content of the p-meson. We obtain in the infrared a 
simple 3 S\ component as a leading component of the p- meson with a 
small admixture of the 3 Z?i partial wave, in agreement with the SU (6) 
flavor-spin symmetry. 



*This paper is dedicated to our colleague Willibald Plessas on occasion of his 60th 
birthday and will appear in special issue of Few Body Systems. 
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1 Introduction 



The SU(6)fs flavor-spin symmetry for the low-lying hadrons in the light 
flavor sector [1] and its roots (as due to the nonrelativistic quark model [2]) 
predated QCD and had numerous phenomenological successes. When QCD 
was established as the fundamental theory of strong interactions it has soon 
become clear that the current quarks of QCD in the u, d sector have tiny 
masses, of the order of a few MeV at the renormalization scale of 1-2 GeV. 
They are very far away from the rather heavy constituent quarks of the 
quark model. In view of this the QCD Lagrangian has approximate chiral 
symmetry. 

At the same time it was clear that this approximate chiral symmetry is 
dynamically broken in the vacuum and this breaking is a source of mass of the 
low-lying hadrons. Different microscopical models (with varying definitions 
of the quark mass) exist for chiral symmetry breaking in the vacuum and they 
indicate that indeed at large space-like momenta the quark mass matches its 
bare values of a few MeV, while at low momenta it runs to the value of the 
constituent quark mass in the definitions given by the models. There also 
exist lattice determinations of such masses at low momenta, see, e.g., [3111], 
though such determinations are manifestly gauge dependent. By fixing a 
gauge and considering a single quark propagator in the background gluonic 
field it is not clear a priori how the confining gluodynamics (that drives 
the dynamics in the color-singlet hadrons) is taken into account. Then it is 
interesting to see how it would be possible to provide a bridge from QCD in 
the infrared to the language of the quark model in a model-independent and 
gauge-invariant manner. 

There is a systematic method to study the hadron composition on the 
lattice - the variational method [5], [6] . In that approach one selects a set of 
interpolating operators {0\, O2, On} with the proper quantum numbers 
that couple to a given hadron. One computes the cross-correlation matrix 

Qjit) = (Oi(t)Oj(O)). (1) 

Masses of the ground and excited states of hadrons with fixed quantum num- 
bers can be extracted from the t-dependence of the eigenvalues of this matrix 
at large Euclidean times t. If the set of operators Oi is complete enough, then 
the eigenvectors represent the "wave function" of the hadron (see the cau- 
tionary remarks at the begin of Sect. H] and the definition (j!9p ). Of course, 
hadrons contain many different Fock components. Our task is to reconstruct 
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the leading one, the quark-antiquark component of the low-lying mesons. 
For this one needs a set of operators that allows one to define uniquely such 
a component. This set of operators must be complete in the qq space with 
regard to the chiral basis. 

All possible qq interpolators for non-exotic mesons in the u, d sector 
have been classified according to transformation properties with respect to 
SU{2) L x SU{2) R and U(1)a chiral groups If one assumes that there is 

no explicit excitation of the gluonic field with the non- vacuum quantum num- 
bers, which is certainly true for low-lying hadrons, then the SU (2) l x SU(2)r 
representations for the quark-antiquark system specify a complete and or- 
thogonal basis. Consequently a set of interpolators that is in one-to-one 
correspondence with all possible chiral representations of SU (2) L x SU (2) R 
is a complete one and can be used to define the qq component of a meson. 
The cross-correlation matrix with such a set can be used to reconstruct the 
qq Fock component. The eigenvectors of this correlation matrix represent 
the qq content in terms of different chiral representations. Observing a su- 
perposition of different chiral representations implies that chiral symmetry 
is broken. 

It turns out that it is also possible to reconstruct a composition of the 
qq component in terms of the 2S+1 Lj basis, where J = L + S are standard 
angular momenta. Indeed, the complete and orthogonal chiral basis can be 
related, through a unitary transformation, to the complete and orthogonal 
2S+1 Lj basis in the center-of-momentum frame [3]. Then diagonalizing the 
cross-correlation matrix with interpolators that span a complete set of chiral 
representations and using this unitary transformation to the 2S+1 Lj basis 
one can obtain from the eigenvectors of the correlation matrix a partial wave 
decomposition of the qq component. 

In QCD the decomposition of a hadron should depend on the scale at 
which we probe this hadron. In other words, what we see in our microscope 
depends on the resolution. In our case "the microscope" is our interpolating 
operator Oj that creates the hadron from the vacuum. A true point-like 
source would correspond to the point-like lattice interpolator in the limit of 
the lattice spacing approaching 0, a — > 0. The point-like interpolator applied 
on the lattice with the spacing a probes the hadron at the scale specified 
by a. In the continuum limit it becomes the true point-like operator and 
probes the hadron at the scale fi 2 — > oo. Here we want to study the hadron 
structure at the infrared scale, where the mass is generated. This scale is 
determined by the hadron size, of the order 0.3 - 1 fm. In lattice simulations 
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we cannot use such a large a, because then the lattice artifacts are too large 
and matching to continuum theory is lost. However, such a low scale can 
be fixed by the gauge-invariant smearing of the interpolators. If we use the 
interpolating operator smeared over the size R in the physical units such 
that R/a ^> 1, then even in the continuum limit a^Owe probe the hadron 
structure at the scale R. Changing the smearing size R we can study the 
hadron content at different scales of the continuum theory at a — > 0. 

In this paper we expand our results on the chiral and partial wave de- 
composition of the p-meson in the infrared, presented in a recent letter [ID] , 
discuss all required details of the formalism and physical interpretation as 
well as give some additional numerical results. 



2 Chiral classification of the quark-antiquark 
interpolators 

The chiral classification of some of the qq interpolators was done in [TT]. A 
complete classification was performed in [TJ E] and is summarized here. 

We consider the two-flavor mesons. All quark-antiquark bilinear opera- 
tors in the chiral limit can be classified according to the representations of 
the SU{2) L x SU(2) R and U(1)a chiral groups. Consider, as example, local 
interpolators of J = mesons, built from quark isodublets q: 

O v (x) = \q{x)r^ 5 q{x) , (2) 

Of (x) = q(x)q{x) , (3) 

O v (x) = iq(x)^ 5 q(x) , (4) 

°a ( x ) = q(x)rq(x) , (5) 

where f denotes the vector of isospin Pauli matrices. The SU {2) L x SU{2) R 
transformations consist of vectorial and axial transformations in the isospin 
space. The axial transformation mixes the currents of opposite parity: 

O t {x) <-> O f0 (x) (6) 

as well as 

O ao (x)^O v (x). (7) 

Hence the currents ([6]) form the basis functions of the (1/2, 1/2) a represen- 
tation of the chiral SU(2)l x SU(2)r while the interpolators ([7]) transform 
as (l/2,l/2) 6 . 
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Table 1: The complete set of qq states (interpolators) classified according to 
the chiral basis. The symbol <-> indicates the states (interpolators) belonging 
to the same representation of SU(2) L x SU(2) R . 



R 


J = 





J = 


1,3,... 


J = 


2,4,... 


(0,0) 






0J++ 


<-> oj— 


oj— 


<-> 0J++ 


(1/2, l/2) a 


1J-+ 


0J++ 


1J+- 


<-> oj— 


1J-+ 


<-> 0J++ 


(1/2,1/2), 


1J++ <-> 


oj-+ 


1J— 


<-> oj+- 


1J++ 


<-> oj-+ 


(0,1)0(1,0) 






1J— 


<-> 1J++ 


1J++ 


<-> 1J— 



As another example consider local interpolators of the p-meson. There 
exist two different bilinear operators with the p-meson quantum numbers 
that have radically different chiral transformation properties. The first one 
is the standard vector current, 

0^(x) = q(x)Yrq(x), (8) 

and the second one is the pseudotensor operator, 

T p (x) = q(x)a°*Tq(x). (9) 

The axial SU(2)l x SU(2)r transformation mixes the vector current with 
the axial vector current, 

O ai (x) = q(x) 1 t l 5 rq(x), (10) 

while the pseudotensor interpolator gets mixed with the interpolator of the 
hi meson, 

O hl (x) =e ijk q(x)cr jk q(x). (11) 
Consequently the operators 

O v p {x)~O ai {x) (12) 

form a basis of the (0, 1) © (1,0) representation of SU(2) L x SU(2) R , while 
the interpolators 

0?(x)~O hl (x) (13) 
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transform as (l/2,l/2) b . A complete set of representations of SU(2) L x 
SU(2)fi and the corresponding quantum numbers are listed in Table [TJ 

In Table [1] the index R determines a representation of the SU(2)l x 
SU(2) R with R = (0,0), (l/2,l/2) a , (l/2,l/2) 6 , or (0,1) © (1,0). All 
states (interpolators) are uniquely specified by the set of quantum numbers 
{R; IJ PC } where we use the standard notations for the isospin /, total spin J 
as well as for the spatial and charge parities PC. The chiral basis {R; IJ PC } 
is obviously consistent with Poincare invariance. The symbol «-> indicates 
that both given states (interpolators) are members of a particular chiral 
multiplet, that is they transform into each other upon SU(2)l x SU(2)r. 

For a particle with the p-meson quantum numbers the (0, 1) © (1, 0) and 
the (1/2,1/2)6 representations are a complete set. Both interpolators, O v p 
and OT, may create the p-meson from the vacuum. If both couple, indeed, 
this signals chiral symmetry breaking. 



3 Transformation from the chiral to the an- 
gular momentum basis 

The chiral basis can be related, through a unitary transformation, to the 
{/; 2 +1 Lj} basis in the center-of-momentum frame [5]: 



\R; u pc ) = E E x^f^cli^cLUi-, 2S+ %> , (M) 



A 9 A- 



where the summation is implied in helicities of the fermion X q and antifermion 
Xq as well as in the orbital angular momenta L such that (— 1) L+1 = P. The 
total spin S is fixed by the quantum numbers IJ PC . Coefficients X\ \- can be 
extracted from Table 2 of Ref . [H] . It follows immediately from Eq. (I14p that 
every state (interpolator) in the chiral basis is a fixed (prescribed by chiral 
symmetry and unitarity) superposition of allowed states in the {/; 2S+l Lj} 
basis. For instance, there are two kinds of the vector states (interpolators) 
with the quantum numbers of the p-meson, which are represented by two 
orthogonal fixed combinations of S- and D-waves: 



|(0, 1) © (1,0); 1 1— ) 
|(l/2, 1/2),; 1 l-> 




U- 3^ (15) 
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with 




In terms of the quark-antiquark bilinears the state |(0, 1) © (1,0); 1 1 ) is 
given by the spatial components of the standard vector current, O v p from (jSJ), 
while the pseudotensor interpolator T p in (JUj) represents | (1/2, 1/2);,; 1 1 ). 
Consequently, diagonalizing the cross-correlation matrix with the O v p and 
T p interpolators and using the unitary transformation we can reconstruct 
the partial wave content of the qq component of the p-meson. 

4 The variational method and the "wave func- 
tion" of a hadron 

The composition of hadronic states in quantum field theory is a subtle issue. 
Whereas in non-relativistic approaches the notion of a wave function and a 
complete basis of states is well-defined, in QFT beyond the ground state there 
is no well-defined single hadron, the state is always a scattering state with 
superposition of many particle components. A given hadron interpolator 
couples in principle to all states with its quantum numbers. 

In the lattice formulations this situation is further complicated by the 
coarseness of the lattice and its symmetries, which mostly have no one-to- 
one correspondence to the continuous space-time symmetries. Lattice in- 
terpolators for mesons, in their simplest disguise, are point like or extended 
color singlet quark-antiquark combinations with suitable combination of their 
Dirac components. 

Combinations of such interpolators may then be formed such as to rep- 
resent lattice analogues of spin representations [12j [T3] . Usually the even 
or odd spin representations mix among them, thus, e.g., a lattice interpola- 
tor representing a J = 1 vector will have contributions which correspond to 
higher odd J states and thus corresponding signals in the correlation func- 
tion. Suitable combinations may improve the situation, but much higher 
statistics is then necessary. This is one of the reasons why lattice studies of 
excited states are progressing slowly. 

1 It is actually possible to construct interpolators with derivatives that also transform 
according to the representations above [8] . 
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The variational method [5J E] (for a recent, more complete set of references 
cf. |14j ) provides a method to find the optimal combinations. Starting with a 
set of interpolators {Ox, O2, ■ ■ ■ , On} one determines the correlation function 
(JTJ) and solves a generalized eigenvalue problem (see below). Given a large 
enough basis set of interpolators the eigenvalues may then be related to the 
eigenstates of the Hamiltonian and the corresponding energy values. 

In order to avoid possible prejudices one should provide a large basis of 
interpolators. On the other hand, using too many such interpolators may 
increase the statistical noise in the analysis. Since the optimal combination 
models the hadron structure it is helpful to be guided by intuition - which is, 
however, another word for prejudice. Thus the selection of a suitable basis 
is close to being an art and one has to balance quality with quantity. 

The normalized physical states \n) propagate in time with 

(n(t)\m(0)) = 5 nm e- Ein)t . (17) 

The interpolating (lattice) operators Oi(t) are projected to vanishing spatial 
momentum and are usually not normalized. 
We compute the correlation function 

00 

C(t)« = (0^)01(0)) = £ a^af >V* w) * , (18) 

n 

with the coefficients giving the overlap of the lattice operator with the phys- 
ical state, 

= (OjO» . (19) 

For interpolating operators Oj spanning an orthogonal basis these values 
would indeed constitute the wave function of state |n) in that basis (more 
exactly, the Bethe-Salpeter amplitude). In this work we are using the term 
"wave function" for this matrix eleemnt.. 

We assume that the correlation matrix ffl8l) can be approximated by a 
finite sum over N states and denote this approximation by 

N 

C(th = £ ai B ^ B) V* (B)t . (20) 

n=l 

The eigenvector and eigenvalues of C will not exactly agree with those of C. 
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It can be shown [5J E] (for a recent discussion see [T51 [33] ) that the gen- 
eralized eigenvalue problem (summation convention) 

C{t) l3 uf ] = X^{t,t )C{t ) l] uf ) (21) 
allows to recover the correct eigensystem approximately, i.e., 

\M(t,t ) = e - E(n) ^ (l + O ( e - A£(n) (*-*°))) . (22) 

Here AE^ may be as small as the distance to the next nearby energy level. 
In [13] it was pointed out that in an interval to < t < 2 1 these contributions 
are suppressed and leading terms even have AE^ n > equal to the distance to 
the first neglected energy level E^ N+1 \ At to ah eigenvalues are 1 and the 
eigenvectors are arbitrary. 

Inserting (|20|) into (l2Tj) one sees that the eigenvectors of the generalized 
eigenvalue problem come out orthogonal to the original wave functions a^ n \ 

N 

(„(»), a M) = £<4 n) X M = c^S nm , (23) 

i=l 

and approximate the correct ones. The normalizing factor we can get 
rid off. 

Defining a sum of lattice operators 

N 

rjW^u^Ot, (24) 

i=l 

we find 

N N 

(m|^)|0> = E«! n) *(0|Ol|m) = ]>>! n) \ (m) = c^5 nm . (25) 
i=i i=i 

Therefore 77 creates a physical state, 

^(n)t|o) = c (»)*| n ) . (26) 

The eigenvector coefficients give the composition of the eigenstate in terms 
of the (non-orthogonal) interpolating operators: 

1 1 N 

at ] = (0\O t \n) = —-(0\O lV ^\0) = -- W^O) . (27) 

c 3=1 

9 



They would agree with Uj if the interpolators were orthogonal, which they 
are usually not. 

The eigenvectors of the generalized eigenvalue problem for hermitian ma- 
trices C(t), (7 (to) obey the orthogonality relation 

(u^\C(t)u^) oc5 nm . (28) 

Indeed, with (1181) and (|23|) we find for large t (summation convention) 



) w (n) 



~ c^*a^e- E(n)t , 



2 e -_B( n )t 



and for the ratio 



(n) 



(») 



(29) 



(30) 



Assuming asymptotically leading exponential behavior this allows to read off 
|dj- | in the asymptotic region. Ideally should not depend on t, in actual 
calculations, however, one identifies a region of t-value where it is compatible 
with a constant. 

We may utilize this result to determine ratios of couplings of the different 
lattice operators to the physical states, 



w&(t) C(t) ijU 



(n) 



» 



(31) 



One computes the ratio for several values of t and identifies its value in a 
plateau region. The ratio tells us how much different interpolating operators 
contribute to the eigenstate |n). This can be used to discuss contributions 
of, e.g., different representations of the vector meson channel, as we do here. 

In a realistic simulation many physical states may contribute and one has 
to use the discussed techniques to single out couplings to the state of interest. 
In some approximation one may instead look just at the ratios of individual 
entries of the correlation matrix, like 



Cu^/C^t) or Caityd^t) 



(32) 
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in the asymptotic t-region, where the excited states contributions are sup- 
pressed pn cm us]. 

The ratio 



(n) 


2 


(0\Oi{t)\n) 


(m) 




(0\Oi(t)\m) 



(33) 



tells us how much the interpolating operator contributes to the eigenstates 
|n) and \m). This can be used to discuss the ratio of decay constants of 
various excitations as done in [T9j [T5 | [T4"]. 



5 Vector meson couplings for local interpola- 
tors 

Here we quote definitions of the local coupling constants of the vector and 
pseudotensor currents to the p-meson in continuum and discuss their relations 
to the matrix elements obtained in the previous section. 

In Minkowski space the corresponding amplitudes are given as 

(0|?(0)yg(0)|V(p;A)) = m,/^, (34) 
<0| (g(OK^O)) (»)\V(p; A)) = ift{p)e^ - e&>) . (35) 

Here V(p; A) is the vector meson state with the mass m p , momentum p 
and polarization A. The vector current is conserved, consequently the vector 
coupling constant is scale- independent. The pseudotensor "current" is not 
conserved and is subject to a nonzero anomalous dimension. Consequently 
the pseudotensor coupling /J(/x) manifestly depends on the scale fi. In the 
rest frame the ratio 

ft (0|7(0)Vg(0)|y(A)) . , 

ft(v) (0\(q(0)^q(0))^)\V(X)) { ' 

coincides with the ratio of matrix elements fl3TT) . 

That ratio can be extracted from the ratio of the vector-vector and 
pseudotensor-vector zero-momentum correlators at large Euclidean times, 
when the excited states do not contribute any more. If the source is located 
at the point (t,x) = (0,0) and the sink is at the arbitrary point (t,x), then 
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at large t one obtains asymptotically 

5>| (q(t,x)Yq(t,x)) (g(0) 7 i g(0)) t |0) ~ (fj) 2 m 2 p exp(-m p t) , (37) 

X:(0|(g(t,f) ( x ^(t,f))( / i)(g(0)7 i g(0)) t |0) ~ /;/J(/i)m>xp(-m p ^8) 

The ratio of these correlators is equal to the ratio fY / 'fT(fj). With the vari- 
ational method the ground state can be identified already at small Euclidean 
time distance and thus the quality of the result can be improved, see Sect. [71 
In [TnidTldB] this ratio has been used to extract fY/fJfa) and to relatdj 
it to the continuum MS-scheme at a scale of 2 GeV. While this method is 
suitable to study the ratio for the ground state and the result agrees with 
that extracted from the variational method (to be compared later on), only 
the variational method can be used to study such ratios for excited states. 

6 Smearing of the interpolators and the res- 
olution scale 

So far we discussed the cross-correlation matrix as obtained with the local 
interpolators. In lattice simulations the hadron interpolators can be built 
with different spatial extent, e.g., using so-called smeared quark sources. For 
example, an isovector interpolator may have the form 

d(x') Sl, T S 0x » u(x") , (39) 

where we omit the Dirac indices; T denotes some Dirac matrix and S xy is 
some gauge transporter from y to x. For the local interpolator at the origin 
one has x' = x" = 0. Summation over x' and x" may be used to define quark 
source smearing, thus improving the signal when computing correlators of 
such interpolators. 

2 Here there is a subtle point. In the MS-scheme of continuum theory one uses a 
renormalization group equation obtained to typically two or three loops to relate the ratio 
at different scales. Such a renormalization group equation cannot adequately represent 
the physics related to chiral symmetry breaking, because chiral symmetry breaking is 
intrinsically a nonperturbative effect. If, at given scale, the p-meson couples strongly to 
both, vector and pseudotensor interpolators, then chiral symmetry is strongly broken in 
the p-meson at that scale. If chiral symmetry is broken, however, it is unclear whether 
perturbative RG is applicable to relate the ratio fY //J(/i) at different scales. 
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One purpose of such a smearing is to improve the quality of the signal 
from the state of interest in order to extract, e.g., its mass. Indeed, with 
the local interpolator the coupling of this interpolator to the physical state is 
determined by the behavior of this state at the origin. With the local inter- 
polators we therefore probe the hadron wave function at the scale fixed by 
the lattice spacing a. If we smear the local interpolator in a gauge-invariant 
way over a spatial region of size R, then the coupling of our interpolator to 
the physical state may be better and the quality of the signal is improved. 
So the smearing typically plays a technical role. In our case, however, we 
give the smearing width a fundamental role - it defines a resolution scale 
at which we study the hadron wave function. Of course, different smearing 
methods may lead to different definitions. 

Only a few quantities in QCD do not depend on the scale. For example, 
the lattice spacing a fixes the ultraviolet cut-off (i.e., the renormalization 
scale) , and the observables such as the ratios of masses of different hadrons, 
their electric charges should not depend on this scale. Most of the quantities 
in QCD do depend on the resolution scale at which they are studied. We want 
to study the hadron wave function in the infrared (where mass is generated), 
i.e., at the very low resolution scale characterized by the typical hadron size. 
Certainly we cannot chose a to be so large since then we lose matching to the 
ultraviolet (continuum) limit of QCD. However, even if we use a reasonably 
small a we can fix a scale where we study the hadron at the smearing size 
R. If R 3> a, it is the size of the smearing R that defines a resolution scale 
where we probe the hadron properties. Consequently, the smearing plays a 
rather fundamental role - it defines a scale at which we study the content 
of our hadron. Physically it means that given a source smearing size R, we 
cannot resolve details of our wave function with the smaller size. Physical 
(continuum) results can be deduced from the extrapolation to the a = 
point while keeping R fixed in physical units. 

We fix the resolution scale in the following way. We substitute a local 
interpolator at the point (t, x) by the interpolator with the same quantum 
numbers but with quark fields smeared in spatial coordinates over the size R 
in a gauge- invariant way around the point (t, x) . The profile of the smeared 
quark fields should be approximately Gaussian with the width R. There- 
fore we use the so-called Jacobi smearing j20l |2T] . A point-like source So is 
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smeared out by acting with a smearing operator M, 



N 



S = MS , M = («#) 



(40) 



n=0 



where if is a hopping term, 
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H = £ [Ufa t) 8 S+M + U](x- 1 t) 4 



(41) 



The smearing extends only over individual time slices, i.e., t is fixed. The 
parameters k (hopping parameter) and N (number of smearing steps) are 
tuned to get an approximately Gaussian shape of the quark source with a 
certain width R in physical units. Different smearing algorithms may be used 
for different definitions of the resolution scale. 

7 Simulation details and results 

In earlier work we have studied so-called chirally improved (CI) fermions 
[221 [23] in the quenched [21 (23 [2S1 US] and the dynamical [27] case for 
two mass-degenerate light quarks. The gauge field action was the tadpole 
improved Liischer-Weisz action [28]. Table [2J gives some information on the 
runs analyzed here, details on the simulation and other observables (e.g., 
hadron masses) can be found in the original papers. The spatial lattice size 
16 corresponds to a physical size close to 2.4 fm, the temporal size is twice 
as large. 

In the quenched simulations the p cannot decay; in our dynamical ensem- 
bles [27J its mass is also above the decay threshold since extra units of relative 
(quantized lattice) momentum are needed for the decay. Below threshold the 
2-pion states will become important and probably modify our results. How- 
ever, we study here only the quark-antiquark contributions to the p lattice 
state. Adding further, more-quark interpolators may affect the overall nor- 
malization of the quark-antiquark component but is unlikely to change the 
ratios of the vector vs. pseudovector operators of the quark-antiquark con- 
tribution. Only this ratio is important for the partial wave decomposition. 

In our analysis we use the variational method, with quark-bilinear meson 
interpolators. We use Jacobi smearing of quark sources with the values for 
k and N such as to obtain a narrow (index n) and a wide (index w) source 
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Table 2: Specification of the data used here; for the gauge coupling only 
the leading value Plw is given, m denotes the bare mass parameter of the 
Cl-action. Further details like the determination of the lattice spacing and 
the 7T- and p-masses are found in [131 12Z]. For the quenched case and for 
the dynamical ensemble A we used 100 configurations, for sets B and C we 
analyzed 200 configurations each. The lattice size is 16 3 x 32 throughout. 



Data 


/3lw 


arriQ 


a [fm] 


77V [Me V] 


m p [MeV] 


Quenched 


7.90 


0.04-0.20 


0.1480(10) 


475-1053 


912-1251 


dyn.: A 


4.70 


-0.050 


0.1507(17) 


526(7) 


922(17) 


dyn.: B 


4.65 


-0.060 


0.1500(12) 


469(4) 


897(13) 


dyn.: C 


4.58 


-0.077 


0.1440(12) 


318(5) 


810(28) 



with effective smearings widths [25] of 0.27 fm and 0.41 fm for the quenched 
ensembles and 0.27 fm and 0.55 fm for the dynamical ones, respectively 

H3H2H. 

The variational method not only allows to identify the excited state (s) - 
depending on the number of interpolators used and the statistical quality of 
the data - but also gives increased stability in the ground state signal. Here 
we discuss only results obtained for the ground state in the isovector vector 
channel, the p-meson (J PC = 1 ). We include the operators 

0\ = u n fd n , 0^ = u w yd w , Op = u P Ydp , 
Ol = SvrV d„ , Ol = u w ^Yd w , Ol = Up^fdp , (42) 

where 7* is one of the spatial Dirac matrices, ^ t is the 7-matrix in (Euclidean) 
time direction, and the subscripts w and n denote the two smearing widths 
of the quark sources, whereas p indicates point quark sources. We denote 
the ratios of the coupling of the p-meson to the different interpolators by 
[a], /at"] (77) for the narrow quark sources and analogously for the other two 
cases w and p. 

In Fig. [1] we compare the ratios (which we call a-ratios henceforth) of 
the coefficients determined from diagonalization of the correlation matrix 
between interpolators O^, O^, O^, as discussed in (13TT) with the direct 
ratio between matrix elements (called d-ratios) according to (1321) . 
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Figure 1: We compare for oX/a^ the a-ratios (left-hand colum) determina- 
tion with the d-ratio (right-hand column) determination for dynamical and 
some quenched data sets. The error bars have been determined with single- 
elimination jack-knife. The horizontal lines indicate the plateau fit range and 
value used in Fig. [31 The data correspond to ensembles given in Table [21 
for the quenched ensemble we show only the results for valence masses 0.04, 
0.08 and 0.16 (in lattice units). 
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The plateaus of the a-ratios are remarkable wide and stable for large 
quark masses and decrease in quality towards smaller masses. However, even 
for run C with m n close to 320 MeV we still observe good quality plateaus 
for time separations At = 3 ... 8. The direct d-ratios have larger errors 
and much worse plateau behavior; for run C only the range At = 5 ... 8 is 
acceptable. The contamination with excited states at smaller distances is 
obvious. The values are, however, compatible with the a-ratios. We also 
confirm that, e.g., the d-ratios {0\0\) j (0^0^) are approximately equal to 

{(oloDKoloi)) 2 . 

The data for the wide sources are qualitatively similar. We therefore use 
the a-ratio averages of the plateau values from 3 to 8 as our value estimate 
for [aX /aT](n) for the narrow sources and [aX /aT](iu) for the wide sources. 

The quark-propagators for point sources were not available but we do 
have correlation matrix entries between smeared source interpolators and 
point sink interpolators. Due to that limitation we cannot use the variational 
method (we do not have the full correlation matrix) but determine the results 
via d-ratios of entries of the correlation matrix giving [al / a J] (p) . 

In the ratios 

{(O v n O v p )/(O v n O T p )) 2 and {{O v w O v p ) / ' (O v w O T p )f (43) 

ideally the effect of the smeared sources should cancel and the ratios should 
agree. In Fig. [2] we compare them and find agreement in the plateau region, 
which, however is shrinking and hardly identifiable for the run C with smallest 
quark mass. It is interesting to note that the wide source ratios appear to 
have less contamination from excited states. We use these for the values in 
Fig. [3] 

Fig. E] exhibits our results for [ajf/aj] for all three situations of inter- 
polator smearing: w, n and p comparing quenched with dynamical data. A 
systematic dependence on the smearing scale of the interpolators is obvious. 
We also find that towards smaller quark masses the quenched results (i.e., for 
the lowest valence masses 0.04 and 0.06 in lattice units) significantly deviate 
from the data with dynamical quarks. 

We observe a clear and systematic dependence of aX jaF p on the smearing 
properties of the hadron interpolators, more precisely: the quark sources 
building the interpolators. As discussed in Sect. [6] the dependence on the 
smearing width (ranging from 0.55 fm down to the unsmeared point scale 
a m 0.15 fm) allows to relate the composition of the p on various infrared 
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Figure 2: We compare (for the dynamical data runs A, B and C) ratios of 
correlation matrix entries (see (f43|) ). which in the plateau region should give 
[Op/aT](p)- The results for narrow (open symbols) and wide (full symbols) 
sources agree in the plateau region, but the plateau range shrinks and is 
hardly justifiable for run C. The error bars have been determined with single- 
elimination jack-knife. The horizontal lines indicate the plateau fit range (for 
the w sources) and the fit value used in Fig. [3j 
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Figure 3: We compare our results for [aY/aT] for the three smearing scales: 
w, n and p (no smearing). The open symbols denote quenched data, the 
full symbols the data with dynamical quarks. The error bars have been 
determined with single-elimination jack-knife. 
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resolution scales. A more systematic study of that scale dependence may 
allow to better compare with effective models. 

Physically the ratio encodes how the chiral symmetry is broken in the 
(quark- ant iquark components of the) p wave function at different scales. The 
(0, 1) © (1,0) (i.e., the vector) and the (|, |)b (i.e., the pseudotensor) repre- 
sentations are a complete set for the qq component of the p-meson. Hence, 
chiral symmetry is broken in the p wave function such that the qq component 
is a superposition of both with a relative ratio shown in Fig. [3j There we see 
clear evidence for the resolution scale dependence of chiral symmetry break- 
ing. Extrapolating to physical quark masses (where we have to assume that 
opening the 2-pion decay channel would not affect the result significantly) 
we expect that the ratio varies from « 1.75 for point interpolators (a resolu- 
tion scale given by the lattice spacing a ~ 0.15 fm) down to ~ 1.25 for the 
interpolators built with smeared, wide quark sources, where the resolution 
scale is given from the smearing width R « 0.55 fm. 

Inverting the unitary transformations (fl5 l) -(fT6 l) we conclude that the qq 
component of the p- meson is varying from predominantly 3 Si wave for the 
large R values with increasing admixture of 3 Di towards smaller R, e.g., 

0.995|l; 3 Si) + 0.096|l; 3 £>i) , 

|i; 3 Si>, 

0.998|l; 3 Si) - 0.059| l; 3 L>i) . (44) 

In Fig. | we show the inverse ratio aj/aj. Towards p ~ 1/R — > oo 
the pseudotensor contribution decouples from the p-meson, fj(p — > oo) = 
0, as follows from renormalization group behavior [16] in the asymptotic 
freedom regime. Hence the inverse ratio should approach for R — > 0, fully 
compatible with Fig. HI At this point the partial wave decomposition is 
determined by the (0, 1) © (1, 0) representation alone. For large R the inverse 
ratio saturates with the 3 Si partial wave strongly dominating over 3 D±, with 
little dependence on the quark (or pion) mass. 



1.75 



777 = 1.25 
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Figure 4: We show the dependence of [a^/a^] on the width of the quark 
sources. For the point sources we use the lattice spacing a as source width 
R. Results for dynamical quarks: A (circles), B (squares), C (diamonds). 

8 Conclusion 

In this paper we have defined and calculated the chiral and partial wave 
content of the qq component of the p-meson at different resolution scales. 
We have used a complete (in terms of chiral symmetry) basis of interpola- 
tors that allows to define the chiral symmetry content of the quark- ant iquark 
component of the p-meson. We have studied this in lattice simulations in the 
quenched limit as well as for configuration obtained with n/ = 2 dynamical, 
mass-degenerate light quarks. We have computed the cross-correlation ma- 
trix of the interpolators and applied the variational method for our analysis. 
The eigenvectors of the cross-correlation matrix supply us then with the di- 
rect information about decomposition of the quark-antiquark component of 
the p-meson in terms of different representations of the chiral group. 

Given such a decomposition we were able, using the unitary transforma- 
tion from the chiral basis to the LS'J-basis, to reconstruct the partial wave 
decomposition of the p-meson at different resolution scales. It turns out that 
at low resolution scales R ~ 0.2 — 0.6 fm the quark-antiquark component of 
the p-meson is a strong mixture of two representations of the chiral group 
(0, 1)©(1, 0) and (1/2, 1/2) & and consequently the chiral symmetry is strongly 
broken at these infrared scales. Only at the deep ultraviolet scale asymp- 
totic freedom requires that the composition of the p-meson is be determined 
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by the (0,1) © (1,0) representation alone. Consequently, at low resolution 
R ~ 0.2 — 0.6 fm the p wave function is predominantly 3 Si wave with a 
tiny admixture of the 3 D\ wave depending on the scale. Only in the deep 
ultraviolet the p is given by the fixed \I = l; 3 + ^Jlji \I = if D x ) 

superposition of the S- and .D-waves with a sizeable contribution of the D- 
wave. This explains successes of the SU(6) flavor-spin-symmetry for the p- 
meson, that explicitly relies on the 3 Si content of the p wave function. Note, 
however, that we observe a resolution scale dependence of the composition, 
whereas the quark model does not. 
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